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Summary. — Apparently, Gamma-Ray Bursts (GRBs) are all but standard can- 
dles. Their emission is collimated into a cone and the received flux depends on the 
cone aperture angle. Fortunately we can derive the aperture angle through an achro- 
matic steepening of the lightcurve of the afterglow, and thus we can measure the 
"true" energetics of the prompt emission. Ghirlanda et al. (2004) found that this 
collimation-corrected energy correlates tightly with the frequency at which most of 
the radiation of the prompt is emitted. Through this correlation we can infer the 
burst energy accurately enough for a cosmological use. Using the best known 15 
GRBs we find very encouraging results that emphasize the cosmological GRB role. 
Probing the universe with high accuracy up to high redshifts, GRBs establish a 
new insight on the cosmic expanding acceleration history and accomplish the role of 
"missing link" between the Cosmic Microwave Background and type la supernovae, 
motivating the most optimistic hopes for what can be obtained from the bursts 
detected by SWIFT. 
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1. Introduction 

The huge emitted power of Gamma Ray Bursts (GRBs) makes them detectable at 
any redshift z < 20. R is therefore natural to use them as probes of the far universe, 
to explore the epoch of reionization (see e.g. 0, 52) an d to study the properties of 
the material between them and us, through the study of absorption lines in the spectra 
of their prompt emission and afterglow. Furthermore, as it seems very likely, GRBs 
are intimately connected with the formation of massive stars, and through GRBs we 
hope to shed light on the star formation history at unprecedented redshifts, without 
the limitation of extinction which may plague samples of distant galaxy resulting from 
deep surveys (0, see also Firmani et al. in these proceedings). However, until now, 
the hopes to use them as standard candles to measure the geometry and kinematics of 
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Fig. 1. - Residuals of the distance moduli of SNe la and of the 15 GRBs with known redshift and 
collimation angle with respect to the case f2M=0.27, Qa=0.73. Also shown are the differences 
of various other cosmological models with respect to the S7m=0.27, J7a=0.73 case: red solid 
line: Qm=0.2, fiA=0.8; orange solid line: fiM=0.37, fiA=0.85; blue dot-dashed line: S7m=0.45, 
Qa=0.95; black dashed line: Qm=1, Qa=0; black dotted line: f2M=0.01, J7a=0.99. The arrow 
marks the redshift of GRB 000131 (2=4.5) which is the most distant GRB of known z. From 
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our universe have been frustrated by the large dispersion of their energetics, even when 
corrected for collimation [3|). This changed with the findings by Ghirlanda et al. 
(2004a, hereafter GGL, see also Ghirlanda et al. in these proceedings) that the spectral 
properties of GRBs are related with the energy radiated within their collimation cones. 
This correlation is tight enough to allow the most optimistic hopes for the cosmological 
use of GRBs in the SWIFT era. This "Ghirlanda" correlation links the energy i? pe ak 
where most of the prompt radiation is emitted with the total energy E-y radiated during 
the prompt phase of a GRB. Since the emission is collimated into a cone, Ej is not 
simply measured by the received fluence. To find it, one needs to measure the semi- 
aperture cone angle, through the achromatic break in the lightcurve of the afterglow. 
The Ghirlanda correlation is very similar to the "stretching" relation of SN la (|13p: 
more powerful SN la have lightcurves that decay more slowly. The Ghirlanda relation 
is then a sort of Cepheid-like relation: even if the collimation-corrected energetics of 
GRBs are not all equal (i.e. they are not standard candles in the strict sense), we have 
a tool to know each of them (GRBs become "known candles"). 

2. The Hubble diagram 

Fig. ^ shows the Hubble diagram in the form of residuals with respect to the specific 
choice of riM=0.27 and f2A=0.73, for SNe la and GRBs, together with different lines 
corresponding to different Qm, Oa pairs. Since it is likely that GRBs follow the star 
formation history, it is also very likely that they exist up to z ^10-20. Even at these 
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redshifts, GRBs are easily detectable. Note that the error bars of GRBs are slightly larger 
than the error bars of SN la, but note also that the maximum redshifts already measured 
of GRBs is z = 4.5, while SN la extends up to a maximum redshift of 1.7. This is also the 
maximum foreseen value of the future SNAP( 1 ) satellite, designed to accurately find and 
observe distant supernovae to accurately constrain cosmological parameters and models. 
Needless to say, it is highly desirable to have more than one class of standard (or known) 
candles, to cancel possible evolutionary effects (which are likely different for different 
classes of objects) and to cure possible extinction effects (for GRBs, the hard X-rays of 
their prompt emission are virtually unaffected by absorption). 

The crucial issue about the use of the Ghirlanda relation for cosmology is what can 
be called the "circularity problem" : to find slope and normalization of the correlation, 
we are obliged to choose a particular cosmology (i.e. a particular pair of CIm, ^a values). 
The correlation has not (yet) been calibrated with low-redshifts GRBs, unaffected by 
the choise of Qm and simply because low redshift GBRs are very few, and we may 
wait for years (even with SWIFT) before building a sample numerous enough of low— z 
GRBs. In the meantime, we have to deal with the circularity problem, and in the next 
section we will describe the 4 existing methods already suggested. 

3. How to cure the circularity problem 

The Ghirlanda relation has already been used to find constraints on the cosmological 
parameters. We briefly describe here the four methods that have been proposed and 
used. In principle, the cosmological parameters to find are not only and ^a, but 
also the parameter w entering in the equation of state of the dark energy P = wpc 2 
where P and p are the pressure and energy density of the dark energy. Note that w 
can be assumed constant or be a function of redshift. One possible parametrization is 
w = wq + w'z. The cosmological A term corresponds to wq = —1 and w' = 0. Consider 
for simplicity this latter case. To the aim to put constraints on Om and Ct\, the basic 
ideas of the four methods are the following: 

1. Dai et al. (2004) assumed that the E pea k-Ey correlation found using a given pair 
of (Om, ^a) values was cosmology independent. In other words, Dai et al. (2004) 
assumed that the slope of the correlation is fixed (they used E 1 oc E^ c ^ k ). Instead, 
this slope corresponds to a specific choice of Qm, ^a- This method could be applied 
if there is a robust theoretical interpretation of the Ghirlanda relation, putting an 
heavy weight to a given slope and normalization of the correlation. 

2. Ghirlanda et al. (2004b, GGLF hereafter), pointing out the circularity problem of 
the Dai et al. approach, proposed instead to use the scatter around the fitting line 
of the Epc^-Ej correlation as an indication of the best cosmology. In other words, 
they find a given Ghirlanda correlation using a pair of 17m, ^a, and measure the 
scatter of the points around this correlation through a y 2 statistics. Then they 
change CIm, ^a, find another correlation, another scatter and another y 2 which 
can be compared with the previous one. Iterating for a grid of CIm, ^a sets, they 
can assign to any point in the CIm, ^a plane a value of x 2 , and therefore draw the 
contours level of probability. 



( 1 ) http://snap.lbl.gov/ 
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3. GGLF also proposed and used the more classical method of the minimum scatter 
in the Hubble diagram as a tool to discriminate among different cosmologies. In 
practice, this method is almost equivalent to the previous one. 

4. Firmani et al. (2005) proposed a more advanced method, based on a Bayesan 
approach. This is the best method to cure the circularity inherent in the use the 
Ghirlanda correlation. The basic novelty of this method is that it incorporates the 
information that the correlation is unique, even if we do not know its slope and 
normalization. The basic idea is to find the Ghirlanda correlation in a given point 
f2 of the (1m, plane, and see how the entire plane (i.e. any other pair of Hm, Oa) 
"responds" to this correlation: the correlation remains fixed, but the data points 
are calculated using the (running) J7m and Rvalues. The scatter of the points 
in the Hubble diagram gives a corresponding \ 2 and we can transform it into a 
probability. Then we calculate the correlation in another Cl point of the plane, and 
repeat the procedure. We repeat that for all (say N) Cl points. At the end we will 
have N sets of probabilities. We can then sum these sets giving a weight to each 
of them. The first time we do that we can assign an equal weight. The sum is a 
"probability surface" characterizing each Cl point. The next time we do this sum, 
we account for this information (i.e. there are certain Cl points more probable than 
others) by assigning a weight to each term of the sum equal to the probability we 
have derived in the previous cycle. Then we iterate until convergence, which is 
reached when the "probability surface" does not change from one iteration to the 
next. An optimization of this method is obtained using a Montecarlo technique. 

To find constraints on wq and to' the approaches are the same as discussed above, but 
bearing in mind that the paucity of data, at present, does not allow to constrain more 
than two parameters, hence in these cases one assumes a flat universe and w' = (and 
derives constraints on J7a and Wq), or the concordance model (to constrain wo and w'). 



4. Cosmological constraints 

In Fig. |21 we show the results taken from Firmani et al. (2005). On the left we show 
the confidence contours on Om, ^a when considering the GRB sample alone (15 objects, 
red lines), SN la alone (156 objects, blue lines) and the combined GRB+SN la sample 
(colored regions). The 68% confidence contours of the 15 GRBs are less restrictive than 
the corresponding SN la contours, but note that the contours of the combined sample 
are now consistent with the concordance model (Qm=0.3, Q\—0.7) within the 68% level, 
whereas the SN la alone are not. The dashed blue lines are the contours corresponding 
to the 14 SN la with z > 0.9. While the extents of the confidence contour regions are 
comparable (the slightly larger errors of GRBs are compensated by their larger redshifts), 
the inclination of these "ellipses" is different, and this is due to the different average 
redshifts of GRBs and this subsample of SN la, as explained below. The right panel of 
Fig. |2Jshows the constraints in the SIm~wo plane, once a flat universe is assumed. Again, 
note that the GRBs alone do not yet compete with SN la, yet the constraints of the 
combined sample makes the concordance model (too = — 1) more consistent, being now 
fully within the 68% contour levels. 
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Fig. 2. - Left: Constraints in the JIm-Oa plane derived from our GRB sample (15 objects, 
red contours), from the "Gold" SN la sample (156 objects; blue solid lines, derived assuming a 
fixed value of Ho = 65 km s" 1 Mpc" 1 , making these contours slightly different from Fig. 8 of 
|15| . and from the subset of SNe la at z > 0.9 (14 objects, blue long dashed lines). The three 
colored ellipses are the confidence regions (orange: 68%; light green: 90%; yellow: 99%) for the 
combined fit of type la SN+GRB samples. Dashed lines correspond to the changing sign of the 
cosmic acceleration [i.e. q(z) — 0] at different redshifts, as labelled. Crosses are the centers of 
the corresponding contours (red: GRBs; blue: SNe la, white: GRB+SN la). The black dot 
marks the ACDM cosmology. The dotted line corresponds to the statefinder parameter r — 1 
(see e.g. |16|h in this case it coincides with the flatness condition. Right: Constraints on wo 
and J1m for a flat cosmology with dark energy whose equation of state is constant (wo = — 1 
corresponds to a cosmological constant). Colored regions: combined (GRB+SN la) constraints 
(color code as in Fig. 2). Dashed lines correspond to the changing sign of the cosmic acceleration 
[i.e. q(z) = 0] at different redshifts, as labelled. Crosses are the centers of the SNe la contours 
(blue: SNe la alone, white: GRB+SN la). The black dot marks the ACDM cosmology. The 
dotted line corresponds to the statefinder r — 1. From |7j. 



5. The Cosmic Whirl 

It is instructive to see the behaviour of the luminosity distance ^(^j^Mj^a) in the 
f^M - plane, for different redshifts. To this aim, we show in Fig. [3]the "stripes" within 
which di, changes by 1% in the f^M - ^A plane, for different redshift, assuming that each 
stripe passes through fiM=0.3 and Ua=0-7- Note the following: 

• The width of the stripes decreases for larger redshifts. This is a consequence of 
the topology of the surfaces of constant g?l : at low redshift the surface is a gently 
tilted plane, at high redshifts the surface is more warped, and there appears a 
"mountain" with the peak close to f^M^ 0.1 and 1. This is shown in Fig. [3| 

for z = 1. 
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Fig. 3. - Left: The Cosmic Whirl. Each stripe represents the loci of points of luminosity 
distance differing at most by 1% to the value calculated for f2M=0.3, f2A=0.7, for different 
values of the redshift, as labelled. Right: The surface of equal luminosity distance calculated 
for 2 = 1. 



• As a consequence, the "stripes" at high redshifts are curved, and at very high 
redshift they surround the "mountain peak" . 

• Note that the width of the stripes at high redshifts become narrower for large 
f2A - values, as a consequence of the increasing slope of the plane. 

This example shows that if we have a sample of standard candles characterized by a 
small range in redshifts, the corresponding confidence contours in the f^M - plane will 
be elongated in the the direction of the stripe of the average redshift of the sample. This 
easily explains why the confidence contours derived by using SN la are elongated in the 
direction of the di,(z ~ 0.5) stripe, while the contours derived by the WMAP satellite 
are elongated along d\,(z ~ 1100) stripe. The contours derived using our GRB sample, 
of larger average redshift than SN la, are more vertically elongated. Therefore there is a 
counterclockwise rotation of the confidence contours when increasing the average redshift 
of the standard candles we are using: only a class of sources spanning a large range of 
redshifts can give an accurate value of Om and £!a- 

6. The future 

The constraints posed by WMAP (|T7|) and SN la [E]), together with cluster 

of galaxies (P) have convincingly led towards the so called "concordance cosmological 
model", a flat universe with fiM=0.73 and f2A=0.27. So, why bother? We care because 
even if we pretend to know how the universe is now, we do not know how it was. This 
ignorance is one of the obstacles towards the understanding of what is the dark energy. 
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Fig. 4. - Contours levels in the !5m, f^Aplane for a simulation of 150 GRBs, according to method 
4 (see text and compare with Fig. 2 showing the results with the current sample of 15 GRBs). 
Blue solid lines are the contours for the 156 SN la of the gold sample of |15| . blue dashed 
lines correspond to the subsample of high-z SN la, colored regions correspond to the combined 
sample of 156 SN la and the 150 GRBs. The blue cross is the position of the minimum y 2 f° r 
SN la only, the red and white crosses are the positions of the minimum for GRB only and 
the combined GRB+SN la sample, respectively. These crosses are partially hidden by the black 
dot, marking the f2M=0.3, Qa=0.7 point. Black dashed lines correspond to the changing sign 
of the cosmic acceleration [i.e. q(z) = 0] at different redshifts, as labelled. 



And this is indeed one of the most important and fascinating challenges of modern (as- 
trophysics. This is why GRBs can be extremely important, because they can bridge 
the gap between relatively close-by supernovae and the fluctuations of the cosmic mi- 
crowave background. They may not be the only tool, since also through the integrated 
Sach-Wolfe effect one hopes to learn how the universe was in the past. But having more 
than one way to measure the universe is certainly not redundant, since each way brings 
its own uncertainties, selection effects and so on. Furthermore, we have ahead a few 
years of SWIFT results, which will hopefully discover a few hundreds of new bursts with 
well measured properties, especially their afterglows at early times, which will greatly 
improve the now poor sample of 15 cosmologically useful GRBs. The fact that the mea- 
sured fluences do not strongly correlate with redshifts (see GGL) means that even for 
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high redshifts GRBs we can well determine the fluence itself and the spectral parameters 
of the prompt emission (i.e. -E p0 ak)- SWIFT itself will follow the early afterglow with 
unprecedented accuracy, and this will certainly help to construct detailed lightcurves (in 
X-rays and in the optical) allowing to pinpoint the jet break time more accurately. For 
bursts at z > 5, therefore in the absence of the optical lightcurve (the optical flux is 
strongly absorbed by Lyman a absorption) one can derive the jet break time from the 
X-ray and the infrared lightcurves, and in this case robotic IR telescopes like REM (for 
the early data, HE|), and larger telescopes (for later data) will help. Fig. ^shows how the 
constraints on the cosmological parameters could improve passing from 15 to 150 bursts, 
a number comparable to the current SN la gold sample of |15j which will be hopefully 
reached after 2-3 years of SWIFT. This simulation has been performed (Lazzati et al. 
in preparation) assuming a "luminosity function" of Ei so and an appropriate redshift 
distribution (see, e.g. H3|). We have assigned relative errors of 20% on ibreak, 16% on 
Speak and 10% on the fluence (these are the average relative errors for the 15 GRBs). 
The average redshift of the simulated sample is (z) ~ 2.26. This is the reason of the 
inclination of the ellipses of the contour levels (cfr. Fig. 0). 

The results shown in Fig. 0] suggest that GRBs, reaching larger redshifts than SN 
la, can constrain £Im, better than SN la. Even tighter constraints can be achieved, 
of course, after a sufficient number of GRBs are observed at small redshifts, allowing 
a cosmology-independent calibration of the Ghirlanda relation, and/or after a robust 
theoretical interpretation is found for this relation. 

For an on-line update: www. merate. mi. astro. it/^ghirla/deep/blink. htm. 

s{c * * 
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